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ABSTRACT 



Context. Chromospheric and prominence plasmas contain neutral atoms, which may change the plasma dynamics through collision 
with ions. Most of the atoms are neutral hydrogen, but a significant amount of neutral helium may also be present in the plasma with 
a particular temperature. Damping of MHD waves due to ion collision with neutral hydrogen is well studied, but the effects of neutral 
helium are largely unknown. 

Aims. We aim to study the effect of neutral helium in the damping of Alfven waves in solar partially ionized plasmas. 
Methods. We consider three-fluid magnetohydrodynamic (MHD) approximation, where one component is electron-proton-singly 
ionized helium and other two components are the neutral hydrogen and neutral helium atoms. We derive the dispersion relation of 
linear Alfven waves in isothermal and homogeneous plasma. Then we solve the dispersion relation and derive the damping rates of 
Alfven waves for different plasma parameters. 

Results. The presence of neutral helium significantly enhances the damping of Alfven waves compared to the damping due to neutral 
hydrogen at certain values of plasma temperature (10000 - 40000 K) and ionization. Damping rates have a peak near the ion-neutral 
collision frequency, but decrease for the higher part of wave spectrum. 

Conclusions. Collision of ions with neutral helium atoms can be of importance for the damping of Alfven waves in chromospheric 
spicules and in prominence-corona transition regions. 

Key words. Sun: atmosphere - Sun: oscillations 
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1. Introduction 

Existence of neutral atoms may change the dynamics of astro- 
physical plasma due to their collision with charged particles. 
Ion-neutral collision may lead to the damping of magnetohydro- 
dynamic (MHD) waves in the solar photosphere/chromosphere 
and prominences (Khodachenko el al. 120041 Forteza et al. 120071 
Zaqarashvili et al. 12011b . Only part of hydrogen atoms are ion- 
ized in the solar photosphere, but the ionization degree increases 
with height due to the rise of plasma temperature. The plasma 
in the transition region/lower corona is almost fully ionized, but 
upper chromosphere, spicules and prominences contain signif- 
icant amount of neutral atoms. Numerous papers study the ef- 
fects of plasma partial ionization in different aspects of the solar 
atmosphere (Haerendel |1992| De Pontieu and Haerendel |1998l 
De Pontieu et al. 120011 Khodachenko and Zaitsev i2002l James 
and ErdelyiEOOll James e t al. l2004l Khodac henko el al. l2004l 
Khodachenko and Rucker 120051 Leake et al. 120051 Leake and 
Arbe r[2006 l Arber et al .[2007l Forteza et al. l2007ll2008l V ranjes 
et al. 120081 Sole r et al I2009al I2009bl Gog oberid ze et al. l2059l 
Carbonell et al l20T0l Singh an d Krish an l20T0l Goodman & 
Kazeminezhad 120101 Tsap et al. 1201 11 1. Most of these papers 
consider only the effects of neutral hydrogen and use a single- 
fluid MHD approximation. Recently, Zaqarashvili et al. (1201 11 1 
showed that the consideration of two-fluid approach significantly 
change the dynamics of MHD waves compared to the single- 
fluid one when wave frequency is comparable or higher than ion- 
neutral collision frequency. They show that the damping rates 



of transverse waves (Alfven, fast magneto-acoustic) due to ion- 
neutral collision reach maximum value at certain wavelength, 
but then begin to decrease for higher frequencies. The two fluid- 
description shed light on few disputed problems, such as, for 
example, the damping of slow magneto-acoustic waves (Forteza 
et al. I2007I I. Therefore, the multi-fluid approach reveals more 
complete dynamics of partially ionized plasmas. 

On the other hand, hydrogen is not the only sort of neu- 
tral atoms, which may be important in the damping of MHD 
waves. The solar plasma may contain significant amount of neu- 
tral helium in certain regions of the solar atmosphere. The up- 
per chromosphere/transition region, spicules and prominences 
are the regions, where the presence of neutral helium may be 
of potential importance for the damping of MHD waves. The 
first step towards the inclusion of neutral helium has been done 
by Soler et al. (120101) . They used the single-fluid MHD approach 
and concluded that the neutral helium causes no significant influ- 
ence on the damping of MHD waves in the prominence plasma. 
However, the effect of neutral helium can be enhanced in the re- 
gions of the solar atmosphere with higher temperature, therefore 
the further analysis is of importance. 

In this paper, we study the Alfven waves in multi-fluid par- 
tially ionized plasma, which contain neutral hydrogen and neu- 
tral helium atoms. The particular attention is paid to the wave 
damping due to ion-neutral collisions and the difference between 
damping rates in partially ionized plasma with and without neu- 
tral helium. We derive the three-fluid MHD equations from ini- 
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tial five-fluid equations and solve the linearized equations in the 
simplest case of a homogeneous plasma. 

2. Main equations 

We study partially ionized plasma, which consists in electrons 
(e), protons (H^), singly ionized helium (He^), neutral hydro- 
gen (H) and neutral helium (He). We suppose that each sort of 
species has Maxwell velocity distribution, therefore they can be 
described as separate fluids. Below we first write the equations 
in five-fluid description and then perform consequent transition 
to three-fluid approach. 

2.1. Multi-fluid equations 

We aim to study the dynamics of Alfven waves, therefore we 
consider the incompressible plasma. We also neglect the viscos- 
ity, the heat flux and the heat production due to collisions be- 
tween particles. Then the fluid equations for each species can be 
written as (Braginskii ll965l Goedbloed & Poedts l2004l l 



V ■ y„ - 0, 



(1) 



«««« \^ + (^«- ^)^«) = -"^p^ - ^««« Ie + ^v„xb\ + 



(2) 



where m^, «„, pa, ¥„ are the mass, the number density, the pres- 
sure, and the velocity of particles a, E is the electric field, B 
is the magnetic field strength, Rg is the change of impulse of 
particles a due to collisions with other sort of particles, e,, - 
+4.8 X 10"'" statcoul is the charge of electrons, protons and 
singly ionized helium (note, that - for neutral particles) 
and c - 2.9979 x 10'" cm s"' is the speed of light. Plasma is 
supposed to be quasi-neutral, which means «e = nH+ + nHe+ ■ The 
description of the system is completed by Maxwell equations 
which have the forms (without displacement current) 



VxE^ , 

c dt 

^ „ 47r . 
VxB^—j, 

c 



where 



(3) 



(4) 



(5) 



is the current density. 

In the case of Maxwell distribution in each sort of particles, 
Ra has the form (Braginskii |1965l l: 



Ra^-Yj^abiVa-Vh), 



where - ata are coefficients of friction between particles a 
and b. 

For time scales longer than ion-electron and ion-ion collision 
times, the electron and ion gases can be considered as a single 
fluid. This significantly simplifies the equations taking into ac- 
count the smallness of electron mass with regards to the masses 
of ion and neutral atoms. Then the five-fluid description can be 
changed by three-fluid description, where one component is the 
charged fluid (electron-Hprotons-Hsingly ionized helium) and the 
other two components are the gases of neutral hydrogen and neu- 
tral helium. 



We use the definition of total density and the total velocity 
of charged fluid as 



P = PH+ +PHe+, 

PH+Vh+ +PHe+VHe+ 



PH+ +PHe+ 



(7) 
(8) 



correspondingly. Then the sum of momentum equations for elec- 
tron, proton and singly ionized helium gives the equation 



dV 1 

P— +P^H+^He+(H'V)H' = -Vp + -jxB+ F,, 

at c 



(9) 



where 



Fi - -(aH+H + ttH+He + ttHe+H + Q'He+He)^+ 



+^H+(Q'He+He + <3'He+H)w - ^He+(a'H+H + 0'H+He)w + 
+ (ffH+H + a'He+H)^^H + (O'H+He + a'He+He)VHe (10) 

is the collision term, p - Pe + Ph* + PHe* is the total pressure of 
charged fluid, = "Ih+^h+Zp, ^He+ = rnue+nue+ /p are relative 
concentrations of ions and w - ^^He+ is the relative velocity 
of protons and helium ions. In derivation of Eq.|9]we neglect the 
electron inertia. Friction of plasma with neutral atoms is mostly 
defined by the ion-neutral collision {ae„/ai„ ~ -\/me/mi <*: 1), 
therefore the collision of electrons with neutrals is also neglected 
here. 

The relative velocity of heavy ions, w is small compared to 
the center of mass velocity. This can be easily seen if one uses 
the momentum equations for proton and singly ionized helium 
without neutral atoms. Subtracting the equations one may get 

dV dw 

(mH+ - mile*)— + (»JH+^He+ + »JHe+^H+)-7- + 

dt dt 



+ ('WH+^He+ - 'niie*^^,)iwV)w 



V/7H+ Vp 



He+ 



«H+ 



+ -WXB- aH+He+ (— + ^— Ih-. (II) 

«He+ / C \«H+ nHe+/ 



The ratio of collision and Lorentz terms is proportional to 
^H+He+ /^H*, where Qh+ = eB/nin+c is the proton gyrofrequency 
(~ 10^ i"' for the magnetic field strength of 10 G) and <5H+He+ is 
the ion-ion collision frequency. This ratio is much smaller than 
unity in the solar chromosphere, therefore the collision term can 
be neglected. The ratio of the second term in the left hand side 
and Lorentz term is proportional to w/Qh+, where w is the wave 
characteristic frequency. This ratio is much less than unity for 
the time scales much longer than ion gyroperiod, therefore the 
inertial term with w can be also neglected. The third term in the 
left hand side is even smaller Then one may get 



(6) IL 



Q.II+ 



(12) 



which yields that w <sc V for the time scales much longer than 
ion gyroperiod. Therefore, the terms in Eqs. l9lfT0l which con- 
tain w, are smaller comparing to other terms. Consideration of 
these terms complicate the calculations, but has no significant 
influence on the final results as we are interested in the effects of 
ion-neutral collisions. Therefore, the terms are neglected in the 
rest of the paper. 

Momentum equations for neutral hydrogen and neutral he- 
lium atoms can be written as 



dVi 



Ph- 



dt 



-Vpii + Fn, 



(13) 
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at 



(14) 



where the collision terms are 

f H - -(ffH+H + ffHe+H + a'HeH)VH + 0'HeHVHe + 
+ (an+H + ffHe+H)!^, (15) 

^He = -(^H+He + ttHe+He + Q'HeH)VHe + aneHVH + 
+ (ttH+He + Q'He+He)^^- (16) 

Ohm's law is obtained from the momentum equation of elec- 
trons and it takes the form 

^/'e 1 . „ 1.7 T> Q'eH++Q'eHe+ . , , 

— + jxB--VxB + — J. (17) 

e«e ce«e c e~n^ 

Here the terms containing the ion relative velocity, w, and 
electron-neutral collision are again neglected. 

Maxwell equation (Eq.|3]l and Ohm's law (Eq.[T7]) lead to the 
induction equation 



dt 



Vx(VxB), 



(18) 



where the battery. Hall and Ohmic diffusion terms are neglected. 



2.2. Ion-neutral collision frequency 

The coefficient of friction between ions and neutrals (in the case 
of same temperature) is calculated as (Braginskii 1 19651 1 



«,n„m,„cr,„- 



(19) 



where T is the plasma temperature, m, (m„) is the ion (neutral 
atom) mass, nii,, = mjm„l{mi + m„) is reduced mass, «,■ («„) is 
the ion (neutral atom) number density, cr,„ = nir, + r„)^ - Anrf 
is ion-neutral collision cross section and k = 1.38 x 10"'^ erg 
K ' is the Boltzmann constant. Using mean atomic cross section 
nrf = 8.7974 x 10"'^ cm^ and the atomic masses of hydrogen 
and helium, the expression can be rewritten for each collision 
pairs as 

an+H ~ 8 X IO^^^zh^ih ^ gcmr^s'\ 

ffH+He ~ 10"-'^«H+«He Vr gcm^-'i"', 
ttHe+H ~ lO^^^nHe+WH Vr g«7r^i"'. 



O'He+He ~ 1.6 X 10 ^^«He+«He ^ g CM 



3 „-l 



(20) 



where T and n are normalised hy I K and cm"^, respectively. 

Collision frequency between ions and neutrals needs addi- 
tional discussion. It is commonly accepted that the collision fre- 
quency of ions with neutrals v,„ = 0^1 111,111 is generally diff'erent 
than the collision frequency of neutrals with ions v„, = ai„lm„n„ 
(De Pontieu et al. 1200 II Vranjes et al. l2008l l. In the first (second) 
case the neutrals (ions) are supposed to have zero initial velocity 
and this approximation leads to the diff'erent values of collision 
frequency. However, mean collision frequency between moving 
ions and neutral atoms is supposed to be the same value due to 
physical grounds. From simple equations of motions of ions and 
neutrals (without pressure and Lorentz terms) one can derive the 
equation for relative velocity between ions and neutrals 



d{Vi - Vn) 

dt 



1 1 

m,n; m„n, 



(Vi - Vn), 



(21) 
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Fig. 1. Height dependence of atmospheric parameters accord- 
ing to FAL93-3 model (Fontenla et al. 119931) . Upper panel: the 
plasma temperature. Middle panel: blue solid line - the ratio of 
neutral hydrogen and electron number densities, «H/«e; green 
dashed line- the ratio of neutral helium and electron number 
densities, «He/«e- Lower panel: the ratio of neutral helium and 
neutral hydrogen number densities, nHe/«H- 



where V, and V,, are velocities of ions and neutral atoms, respec- 
tively. This equation gives the collision frequency between ions 
and neutrals as 



1 



1 



(22) 
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This expression significantly changes the value of collision fre- 
quency between ions and neutrals. For example, calculation of 
collision frequency in the parameters used in Zaqarashvili et al. 
(I20m Eq. 30) gives 16 s^' instead of 4 s^'. 

The chromospheric temperature of 1.6 x 10"* K yields the 
ion-neutral collision frequencies (as expressed by Eq.l22li as: 



VH+H - «H+H 



1 



1 



mil* nil* 



mntiH 



61.4 s- 



1 1 

VH+He - ttH+He I + 



1 1 
VHe+H - ffHe+H I ^ 



19 s-\ 



7.2 s-\ 



^He+He - CHe+He 



1 



1 



+ 



2.8 s- 



(23) 



where we used proton, neutral hydrogen, singly ionized helium 

l"'cm-^ 1.05 X 
(Fontenla et al. 



and neutral helium number densities of 9.09 x 10'" cm ^, 1.05 x 



,10 ,^-3^ g 9 ^ ^^-3 ^jjjj 2 5 X 10'^ cm-3 



10- c m 

[19931 model FAL93-3) 

One may expect that the effects of ion-neutral collision may 
be enhanced for the waves with frequencies near ion-neutral col- 
lision frequency (Zaqarashvili et al. l201 lb . Therefore, higher fre- 
quency waves are probably more affected by ion-neutral colh- 
sions. 

2.3. Partial ionization in the solar chromosphere 

Before the study of damping of linear Alfven waves due to ion- 
neutral collision, let us briefly discuss the plasma ionization in 
the solar lower atmosphere. The plasma is weakly ionized in the 
photosphere/lower chromosphere, but becomes more and more 
ionized with height. The increase of temperature with height 
leads to the ionization of hydrogen and helium atoms. Hydrogen 
and helium are almost fully ionized in the solar corona, so the 
transition occurs near the region of sharp temperature rise. 

Figure 1 shows the dependence of plasma parameters on 
height according to FAL93-3 model (Fontenla et al. ll993l l. This 
model includes the dependence of ionization degree on heights 
for both, hydrogen and helium. The upper panel shows the 
plasma temperature vs height according to this model. The tem- 
perature minimum is located near 500 km above the basic of 
the photosphere, while the transition region is just above 2000 
km height. The neutral hydrogen number density is much higher 
than the electron number density at the lower heights, but be- 
comes comparable near ~ 1900 km, which corresponds to the 
temperature of 9400 K (blue solid line on the middle panel). 
Hydrogen atoms quickly become ionized above this height. The 
neutral helium number density is also higher than the electron 
number density at the lower heights (green dashed line on the 
middle panel). They become comparable near ~ 1600 km, which 
corresponds to the temperature of ~ 7300 K. The ratio becomes 
smaller and smaller just above this height. The important pa- 
rameter for this study is the ratio of neutral helium and neutral 
hydrogen number densities (the lower panel). This ratio stays 
nearly constant (~ 0.1) up to 1500 km height, then it quickly in- 
creases up to 0.2 at height of 2000 km. Therefore, the effect of 
neutral helium in the damping of Alfven waves should become 
important at higher altitudes. 

In the next section we study the damping of linear Alfven 
waves due to ion collision with neutral hydrogen and helium 
atoms. 



3. Damping of linear Alfven waves 

We consider the simplest case of a static and homogeneous 
plasma with homogeneous unperturbed magnetic field, B-, di- 
rected along the z axis. We assume the wave propagation in xz 
plane i.e. d/dy = 0. Then the Alfven waves are governed by the 
equations 



5Vv 



B- dby 



dt 4npo dz 



V,. + Vhv + VHe 

PO Po PO 



dv] 



Hey 



dt 



pm 



^He ^ 
PHeO 



an 



• Q'HeH Q'HeH 

-VHy H VHey, 



Pho 



ttHe + Q'HeH 



Pm 



PHeO 



-^Hey 



Q'HeH , 
PHeO 



-VHy 



db. 



dv, 
B — . 

dt ~ dz 



(24) 
(25) 
(26) 
(27) 



where v,. (vHy, VHey) are the perturbations of ion (neutral hydro- 
gen, neutral helium) velocity, is the perturbation of the mag- 
netic field, and po {pm, PHeo) are unperturbed ion (neutral hy- 
drogen, neutral helium) density, respectively. Here we use the 
definitions 



Q'H - ffH+H + Q'He+H, Q'He - Q'H+He + Q'He+He- 



(28) 



Collision between neutral hydrogen and neutral helium 
atoms, expressed by Q'HeH, does not contribute significantly 
in the damping of Alfven waves. Moreover, for example, 
Q'HeH /Q'HeH* ~ «h/«h+ < 1 in spiculcs and in prominence 
plasma, therefore we neglect the corresponding terms for sim- 
plicity. Then, Fourier analysis with exp[i{k,z - cot)] gives the 
dispersion relation 

^H^HeflHflHero-'* + '[^HflH(l + ^He) + ^HeflHe(l + ^h)]^^^''- 
-[^He + ^H + 1 + ^H^HeflHflHelTO-^ - '[^HAH + ^HeOlie]'^ + 1 = 



= 0, 

where 



w k.VAPo 
an = 



an 



(29) 



kzVAPO . PH . PHe 

flHe , gn - — , gHe - — , 

Q'He Po Po 



Va = 



^|4npQ 



(30) 



The dispersion relation has four different roots: the two com- 
plex solutions, which correspond to Alfven waves damped by 
ion-neutral collision and two purely imaginary solutions, which 
correspond to damped vortex solutions of neutral hydrogen and 
neutral helium fluids. 

In order to study the effects of neutral helium, the damping 
rates of Alfven waves obtained from Eq. (|29T l should be com- 
pared with the damping rates obtained in the absence of neutral 
helium. The dispersion relation corresponding to Alfven waves 
in two-fluid plasma with only neutral hydrogen can be easily ob- 
tained from Eq. (l29l l (see also Zaqarashvili et al. 201 1) 



^nanm + i{\ + ^n)Tn - ^nanm - / = 0. 



(31) 



This dispersion relation has three different roots: the two com- 
plex solutions, which correspond to damped Alfven waves and 
a purely imaginary solution. The purely imaginary solution is 
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Fig. 2. Damping rate (imaginary part of frequency) of Alfven 
waves, m - oj/Lva, vs normalized Alfven frequency, a - 
kzVA/yH, for three deferent values of chromospheric tempera- 
ture. Blue crosses correspond to the damping rates due to ion 
collision with neutral hydrogen atoms only, while red asterisks 
correspond to the damping rates due to ion collision with both, 
neutral hydrogen and neutral helium atoms. Red (blue) solid line 
corresponds to the damping rate derived in the single-fluid ap- 
proach (Eq.[32b with (without) neutral helium. 



associated with the neutral fluid. This is clearly seen if one ne- 
glects the ion-neutral collision and consequently the neutral fluid 
behaves independently. Then, the fluid has vortex, co = 0, so- 
lution in the incompressible limit, which gains imaginary part 
when collision is included. Therefore, the purely imaginary so- 



lution of dispersion relation Eq. OTT i corresponds to damped vor- 
tex solution of neutral hydrogen fluid (Zaqarashvili et al. 201 1). 
Then, it is straightforward that the fourth imaginary solution in 
Eq. ( l29b arises owing to the vortex solution of neutral helium 
fluid as mentioned above. 

Fig. 2 displays the damping rates of Alfven waves in par- 
tially ionized isothermal plasma with (red asterisks, derived from 
Eq. l29b and without (blue crosses, derived from Eq. [STT i neu- 
tral helium vs normalised Alfven frequency for three different 
temperatures in the chromosphere. The Alfven frequency k,VA 
is normalised on the collision frequency between neutral hydro- 
gen and charged heavy particles, vh = vh+h + VHe+H- The values 
of temperature and number densities are taken from FAL93-3 
model (Fontenla et al. 11993b . The lower panel shows the damp- 
ing rates for the temperature of 8100 K with proton, neutral hy- 
drogen and neutral helium number densities of 1 .33 x 10" cm"^, 
3.71 X 10" cm"^ and 4.8 x 10'° cm"\ respectively. The middle 
panel shows the damping rates for the temperature of 10 000 K 
with proton, neutral hydrogen and neutral helium number den- 
sities of 1.07 X 10" cm-^ 8.47 x 10"' cm^^ and 1.13 x 10'" 
cm"^, respectively. The upper panel shows the damping rates for 
the temperature of 16 000 K with proton, neutral hydrogen and 
neutral helium number densities of 9.09 x 10'° cm"\ 1 .05 x 10"' 
cm"^ and 2.5 x lO' cm"^, respectively. The three different tem- 
peratures correspond to 1790 km, 1995 km and 2015 km heights 
from the photosphere. 

According to Fig. 2 the influence of neutral helium atoms 
on the damping of Alfven waves significantly depends on the 
plasma temperature. The effect of neutral helium is less pro- 
nounced for the temperature of ~ 8000 K (lower panel). But 
it gradually increases for higher temperatures (middle and up- 
per panels). The dependence of the damping rate on frequency 
is similar to that for the neutral hydrogen: the damping rate 
increases at lower frequencies, reaches its maximum at partic- 
ular frequency and then begins to decrease for higher values 
(Zaqarashvili et al. 1201 11 1. The damping rate of Alfven waves 
has its maximum near ion-neutral collision frequency. 

The results show that the neutral helium atoms significantly 
enhance the damping of Alfven waves in the chromospheric 
plasma. For example, the upper panel of Fig. 2 corresponds to 
the parameters of spicules with the temperature of 16000 K 
(Beckers |19681 l. The presence of neutral helium enhances the 
damping rate almost three times for low-frequency waves and 
almost twice for high-frequency ones. Therefore, neutral he- 
lium atoms should be taken into account, alongside with neu- 
tral hydrogen, in the study of transverse waves in solar spicules 
(Zaqarashvili and Erdelvi l2009l l. 

Recently, Soler et al. (12010b have studied the damping of 
MHD waves due to ion-neutral collision in the single-fluid de- 
scription of prominence plasma including neutral helium atoms. 
They concluded that the neutral helium has no significant influ- 
ence on the damping of MHD waves. We check the result in 
the three-fluid MHD description taking plasma conditions typi- 
cal for prominence cores with 8000 K temperature. Proton and 
neutral hydrogen number densities are taken as 9 x 10^ cm"^ 
and 2.1 x lO'" cm"^, respectively (Labrosse et al. 12010b . The 
number density of neutral helium is taken as 10 % of the neu- 
tral hydrogen i.e. 2.1 x 10^ cm"^ and the number density of 
singly ionized helium is taken as 2.1 x 10** cm"^. Fig. 3 dis- 
plays the damping rates of Alfven waves in partially ionized 
isothermal plasma with (red asterisks) and without (blue crosses) 
neutral helium. We see that the neutral helium has insignifi- 
cant effect on the wave damping in considering parameters, es- 
pecially for low-frequency waves, as it was shown by Soler et 
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Fig. 3. Damping rate (imaginary part of frequency) of Alfven 
waves, vj - oj/k~VA, vs normalized Alfven frequency, a - 
k~\'A/vu, for the prominence plasma. Red asterisks (blue crosses) 
correspond to the solutions with (without) neutral helium. Red 
(blue) solid line corresponds to the damping rate derived in the 
single-fluid approach (Eq.[32li with (without) neutral helium. 

al. (1201 Oi l. However, the neutral helium can be of importance in 
prominence-to-corona transition regions, where the plasma tem- 
perature is higher than in cores. Solid lines on this figure indi- 
cate the damping rates in the single-fluid approach. The damp- 
ing rates derived in single- and multi-fluid approaches are sim- 
ilar for low-frequency waves, but have completely different be- 
haviour when wave frequency approaches ion-neutral collision 
frequency. 

4. Discussion 

Solar chromosphere/photosphere and prominences contain sig- 
nificant amount of neutral atoms, which may lead to the damp- 
ing of MHD waves. In some cases, the damping can be of im- 
portance and should be taken into consideration. The effect of 
neutral hydrogen on the damping of MHD waves is well stud- 
ied in both, single fluid (Khodachenko el al. 12004! Forteza et al. 
12007b and multi-fluid approximation (Zaqarashvili et al. 1201 II) . 
However, the collision of ions with neutral helium atoms for 
some cases is important as shown here. Hydrogen atoms quickly 
become ionized with increasing temperature, while the num- 
ber of neutral helium may be still significant. This is clearly 
seen in Fig. 4, which is plotted according to the atmospheric 
model FAL93-3. The ratio of neutral helium to neutral hydrogen 
number densities grows significantly in the temperature interval 
10000-40000 K. 

We study the damping of Alfven waves due to collision of 
ions with neutral hydrogen and neutral helium atoms. The three- 
fluid description of plasma is used, where one component is 
electron-proton-singly ionized helium and the other two com- 
ponents are neutral hydrogen and neutral helium gases. The 
dynamics of linear Alfven waves in homogeneous isothermal 
plasma with constant unperturbed magnetic field is considered. 
The exact dispersion relation for the Alfven waves is derived 
and solved for different plasma parameters. The damping rates 
due to ion collision with neutral hydrogen and neutral helium 
atoms are derived and compared with those obtained only for 
neutral hydrogen. The analysis shows that the collision between 
ions and neutral helium can be of importance at certain values 



Fig. 4. The ratio of neutral helium and neutral hydrogen number 
densities, «He/«H vs plasma temperature according to FAL93-3 
model (Fontenla et al. ll993l ). 

of plasma temperature, when the hydrogen already begins ion- 
ization, but the neutral helium is still presented. This happens 
for T ~ 10000 - 40000 K, which corresponds to upper chromo- 
sphere, spicules and prominence-corona transition regions. The 
presence of neutral helium enhances the wave damping by few 
times compared to the damping only due to ion-neutral hydrogen 
collisions (Fig. 2, middle and upper panels). 

The collisions between ions and neutral helium atoms can be 
of importance in solar spicules which correspond to the upper 
panel on Fig. 2. The plasma temperature, electron density and 
ionization degree are almost constant along spicules (especially 
in their upper part), therefore the isothermal plasma is a good 
approximation there (Beckers 119681) . The neutral helium atoms 
may significantly enhance the damping of transverse waves, 
which are continuously observed in spicules (Kukhianidze et al. 
2()()6. De Pontieu et al. 2007 , Zaqarashvili et al. 120071 see also 
recent review of Zaqarashvili and Erdelvi 120091 ). 

The presence of neutral helium has no significant influence 
on Alfven wave damping in prominence cores (Fig. 3). This 
makes still valuable the results of Soler et al. (120101 ). who studied 
the damping of MHD waves in single-fluid approach. However, 
the neutral helium can be of importance in the prominence- 
corona transition region, where the plasma temperature is higher 
than in the prominence cores. 

Damping rates derived in the three-fluid approach reach a 
peak near the ion-neutral collision frequency and then decrease 
for higher frequencies unlike to the single-fluid approach, where 
the damping linearly increases with increasing frequency (solid 
lines on Figs. 2-3, see also Zaqarashvili et al. l201 lb . This may re- 
quire the reconsideration of high-frequency Alfven wave damp- 
ing in the chromosphere. More sophisticated simulations taking 
into account the gradients of temperature and ionization degrees 
with height in the different parts of solar atmosphere are needed 
to be performed. 

In order to understand the wave damping due to ion-neutral 
collision, we consider the simplest case of partially ionized 
plasma consisting of electrons, protons and neutral hydrogen. 
Then, the well-known expression of damping rate in the single- 
fluid approach (Braginskii il965i Khodachenko el al. |2004I 
Forteza et al. 12007b can be rewritten as 

^i^-, = ^ = 777^ (32) 
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where ptot is the total density (ion+neutral) and v,„ is the ion- 
neutral collision frequency defined by Eq. (l22T i. ^„,^, are nor- 
malised by Ptot- This expression clearly indicates that the nor- 
malised damping rate depends on the ratio of Alfven and ion- 
neutral collision frequencies and the ratio of neutral and ion fluid 
densities. Early statement that the damping rate depends on the 
magnetic field strength can be translated as follows: the increase 
of the magnetic field leads to the increase of the Alfven speed, 
therefore the waves with the same wavenumber have higher fre- 
quency, which is closer to ion-neutral collision frequency, hence 
the normalised damping rate increases. In fact, it is the ratio 
of neutral and ion fluid densities, which determines the damp- 
ing rate of particular wave harmonic. High-frequency waves 
are damped quickly. However, this statement is valid for low- 
frequency wave spectrum below ion-neutral collision frequency. 

FAL93-3 model does not include the coronal values of 
plasma parameters and the ionization degree of helium in the 
corona is unknown. However, it is clear that the coronal plasma 
contains only a small number of neutral helium, therefore the ef- 
fects of neutral helium should be smaller in the corona. However, 
certain influence of neutral helium on wave damping in the 
corona is not ruled out, and further analysis is required to make 
a firm conclusion. 
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5. Conclusions 

Existence of neutral helium atoms, alongside with neutral hydro- 
gen, significantly enhances the damping of Alfven waves at cer- 
tain interval of plasma temperature due to ion-neutral collision 
as compared to the case of only neutral hydrogen. Neutral he- 
lium causes the increase of damping rates few times in spicules 
and prominence-corona transition region at the temperature T ~ 
10000-40000 K. 

Alfven waves have maximal damping rates at some fre- 
quency interval peaking near the ion-neutral collision frequency. 
At the same time, the damping rate is reduced for higher fre- 
quencies. 

The study of high-frequency Alfven wave damping in the 
solar chromosphere with realistic height profile of ionization de- 
gree needs to be revised in the future including both, neutral hy- 
drogen and neutral helium atoms. 

Acknowledgements. The work was supported by the Austrian Fonds zur 
Forderang der wissenschaftlichen Forschung (project P21 197-N16). T.V.Z. also 
acknowledges financial support from the Georgian National Science Foundation 
(under grant GNSF/ST09/4-310). The authors thank to the referee for construc- 
tive comments. 



References 

Arber T.D., Haynes, M. & Leake, J.E., 2007, ApJ, 666, 541 
Beckers, J.M., 1968, Solar Phys., 3, 367 
Braginskii, S.L, 1965. Reviews of Plasma Physics, 1, 205 
Carbonell, M.. Forteza, P, Oliver, R. & Ballester, J.L., 2010, A&A, 515, 80 
De Pontieu, B. & Haerendel, G., 1998, A&A, 338, 729 
De Pontieu, B., Martens, PC.H., & Hudson, H.S., 2001, ApJ, 558, 859 
De Pontieu, B. et al., 2007, Science, 318, 1574 
Fontenla, J.M., Avrett, E.H., & Loeser, R., 1993, ApJ, 355, 700 
Forteza, R, Oliver, R., Ballester, J.L. and Khodachenko, M.L., 2007, A&A, 461, 
731 

Forteza, P, Oliver R. & Ballester, J.L., 2008, A&A, 492, 223 

Goedbloed, H. & Poedts, S. 2004, Principles of Magnetohydrodynamics, 

Cambridge University Press 
Gogoberidze, G., Voitenko, Y, Poedts, S. & Goossens, M., 2009, ApJ, 706, 12 
Goodman, M.L. & Kazeminezhad, E, 2010, ApJ, 708, 268 
Haerendel G., 1992, Nature, 360, 241 
James, S.P & Erdelyi, R., 2002, A&A, 393, Lll 
James, S.P, Erdelyi, R. & De Pontieu, B., 2004, A&A, 406, 715 



